1 



Emission angle dependent HBT at RHIC and beyond 
Peter F. Kolb a and Ulrich Heinz b 

department of Physics and Astronomy, SUNY at Stony Brook, NY 11794-3800, USA 
b Physics Department, The Ohio State University, Columbus, OH 43210, USA 

We study the geometrical features of non-central heavy ion collisions throughout their 
dynamical evolution from equilibration to thermal freeze-out within a hydrodynamic pic- 
ture. We discuss resulting observables, in particular the emission angle dependence of the 
HBT radii and the relation of these oscillations to the geometry at the final stage. 

1. Introduction 

The systematic investigation of observables generated by the broken azimuthal symme- 
try of non-central nuclear collisions |1| has lead to major insight in the dynamics of the 
expansion. In particular the large signal of anisotropic particle flow requires a surpris- 
ingly efficient microscopic rescattering mechanisms || and even agrees with the limit of 
infinitely short rescattering lengths, namely hydrodynamic calculations thus pro- 

viding access to the nuclear equation of state under extreme conditions. Furthermore, 
the large signal is critically dependent on rapid thermalization and early pressure in the 
system, at timescales smaller than 1 fm/c @. Experimentally and theoretically even more 
challenging is the extraction and interpretation of the azimuthal dependence of coordi- 
nate space observables such as two-particle correlations [[SHTlH • Its understanding however 
offers an additional handle to study the space-time geometry of the source, helps to un- 
derstand and clarify persisting problems |12| , and ultimately is one more piece in our 



effort to understand the complete dynamics of heavy ion collision. 



2. Fireball initialization and evolution 

To elucidate the connection of these observables and the geometry that they reflect, 
we study two scenarios in the following. In the first one (referred to as 'RHIC1') we 
choose the initial and freeze-out conditions such that the particle yield as well as the pion 
and antiproton spectral slopes agree well with experimental data from central Au+Au 
collisions at 130 AGeV (for details, see The maximum temperature in the center 
of the collision region is found to be To = 340 MeV at an assumed equilibration time of 
t = 0.6 fm/c. The spectral slopes require a decoupling temperature of T dec = 130 MeV, 
where fluid elements are assumed to liberate their particle content boosted by the local 
flow velocity. Geometry, without involving further parameters, then determines the ini- 
tialization of non- central collisions. Experimental results on elliptic flow confirmed the 
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Figure 1. Time evolution of the spatial ec- 
centricity e x and momentum anisotropy e p for 
RHIC1 (solid) and IPES (dashed). 
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Figure 2. Freeze-out hypersurface of the 
two systems. Shown are cuts at y — 
(solid) and x = (dashed). 



large values expected from hydrodynamic predictions and justified this idealized ap- 
proach as long as the produced fireball is sufficiently large to allow for thermalization 
and 'macroscopic behavior', a requirement which expectedly breaks down above a certain 
impact parameter. For Au+Au collisions at ^/snn — 130 GeV the data indicate that this 
limit is reached for b ~ 7 — 8 fm || . 

With the second scenario we want to stress the geometrical features of the source and 
the imprints on experimental observables. To enhance the relevant signals, we employ 
more extreme conditions of To = 2 GeV at tq = 0.1 fm/c and Td oc = 100 MeV. As the 
system stays above the QCD transition temperature for most of the time of the evolution 
we employ an ideal gas equation of state for simplicity (for reasons that will become clear 
later we label this system 'IPES'). At the moment it is not clear whether these conditions 
are accessible at future colliders, but the calculation allows us to trace back the origin of 
qualitative changes in the observables which might already occur at LHC energies. In the 
following we study the evolution of both systems for an impact parameter of 7 fm on the 
basis of an ideal hydrodynamic evolution with longitudinal boost-invariance 0. 

Upon impact the nuclei produce secondaries and large amounts of entropy in their 
geometrical overlap region. This is generally modeled by contributions from wounded nu- 
cleons and/or binary collisions to the initial energy- or entropy density fields in the overlap 
region M. We characterize the spatial eccentricity of the fields by averaging over the trans- 
verse plane (at z = 0) with the energy density as weight. Defining e x = (y 2 — x 2 ) / (y 2 + x 2 ), 
we can follow its evolution through time (see Fig. p. If the impact parameter points in 
x-direction, e x is positive during the initial stages when the overlap-ellipsoid is elongated 
out of the reaction plane. As time evolves, the stronger pressure gradients in the short 
direction force a stronger expansion in plane than out of plane, thereby reducing the ini- 
tial anisotropy and eventually leading to an overshoot such that the source appears in 
plane elongated at the very end of its existence (e x < 0). This is in particular true for 
the 'IPES' case, which exhibits much stronger pressure gradients and therefore stronger 
transverse expansion ('IPES' stands for 'In Plane Elongated Source'). The generation 
of the accompanying anisotropic flow is characterized through the momentum anisotropy 
e p = (T xx — T yv ) I (T xx + T yy ) where the T n are the diagonal elements of the energy mo- 



3 



mentum tensor of the fluid. Ultimately, this momentum anisotropy is reflected as elliptic 
flow. Fig. [I] shows clearly how the geometrical eccentricity causes a rapid buildup of 
momentum anisotropy (the early saturation of e p at about 3 fm/c in the RHIC1 case is 
caused by the mixed phase of the employed equation of state [0]). 

As particle emission occurs throughout the evolution from the surface of the fireball, 
it is not clear what average geometry will be imprinted on the particle distributions. 
Particles with high transverse momenta preferentially emerge from regions with the largest 
collective transverse velocities, that is from the freeze-out hypersurface at its largest radial 
extension, whereas the low p± particles are emitted when the interior of the fireball freezes 
out in a timelike manner (see the freeze-out hypersurfaces in Fig. |^). Therefore particles 
of different p± (or particle pairs with different mean transverse momentum K± as they are 
studied in HBT) probe the system at different times and different regions of the fireball. 



3. Azimuthally sensitive HBT-analysis 

From the hydrodynamically determined freeze-out hypersurface £ we evaluate the den- 
sity of emitted pions of momentum K through S(x, K) oc / s K ■ d 3 a(x') fBE(E(x'),T dec ) x 
<5 4 (x — x') with E(x') = K ■ u(x') where u is the local flow velocity. With this source 
function S(x, K) as a weight we can calculate the average emission points and variances 
and collect them in the spatial correlation tensor S^ u = {x^xj) — (x^)(x u ). Rotating this 
correlation tensor into the direction of the particle pair momentum (angle $), one obtains 
the 'width' and 'depth' of the emission region with respect to the direction of observation 
and can make the connection |J to the experimentally deduced HBT radii M. 

The 'sideward' radius receives contributions only from the (x, y)-components of the 
correlation tensor and thus can be related directly to the geometry of the fireball. More 
specifically it reflects the 'width' of the emission region of particles of similar momentum 
orthogonal to their average momentum (the direction of observation). For large transverse 
momenta, these 'homogeneity regions' are narrow slivers close to the rim of the freeze-out 
hypersurface whose shape is quite different from that of the total source |1C|. The 'out' and 



'out-side' terms mix geometrical and temporal emission information, with the temporal 
part gaining weight with increasing transverse momentum. In the calculations we can 
investigate both parts separately, but experimentally they appear always in combination. 

Fig. [| displays the oscillations of the squared radii for the RHIC1 calculations. Each 
radius shows its own oscillation pattern which remains qualitatively unchanged for all 
transverse momenta. The sideward oscillation reflects the out-of-plane extended source: 
observing such a source from the x-direction it appears to be larger sideways than when 
observed from the y-direction. This remains true even at large K± where the emission is 



concentrated along a thin sliver near the rim of the fireball JED] . The outward oscillations 
confirm this behavior from an 'orthogonal' viewpoint. The sign of the out-side oscillation 
can be traced back to a 'tilt' of the emission region relative to the emission direction 
P!D[ which is generated by the out-of-plane deformation of the source. Even though 
hydrodynamic calculations fail in the quantitative description of the freeze-out geometry 



121 , the overall source shape appears to be in agreement with preliminary data ||. 

Moving to the IPES calculation a drastic change occurs (Fig. ||): The oscillations flip 
sign at K± ~ 0.15 GeV, showing the same oscillation pattern as at RHIC1 for larger 
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Figure 3. Oscillations of the sideward, Figure 4. As Fig. ^ for the IPES system, 
outward, out-side and longitudinal radius The purely geometrical contributions are 
for the RHIC1 system. marked by circles. 



K± but opposite oscillations below. The geometric contributions to the out and out-side 
radius (marked by circles) always oscillate opposite to the RHIC1 case, reflecting the in- 
plane elongation of this source. At larger Kx this naively expected behavior is masked 
by growing temporal contributions, leading eventually to the same qualitative oscillation 
pattern as seen for RHIC4. Because the homogeneity regions do not trace out the interior 
of the source but concentrate near the rim of the system, the oscillations of the sideward 
radius also change sign at larger K±. The true geometry of the source is therefore only 
reflected in the correlations of the particles with the smallest transverse momenta. 
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